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In Poulin and Serpico [Phys. Rev. Lett. 114, 091101 (2015)] we have recently argued that when 
the energy of a photon injected in the primordial plasma falls below the pair production threshold, 
the universality of the nonthermal photon spectrum from the standard theory of electromagnetic 
cascades onto a photon background breaks down. We showed that this could reopen or widen 
the parameter space for an exotic solution to the “lithium problem.” Here we discuss another 
application, namely the impact that this has on nonthermal big bang nucleosynthesis constraints 
from 4 He, 3 He and 2 H, using the parametric example of monochromatic photon injection of different 
energies. Typically, we find tighter bounds than those existing in the literature, up to more than 1 
order of magnitude. As a consequence of the nonuniversality of the spectrum, the energy dependence 
of the photodissociation cross sections is important. We also compare the constraints obtained with 
current level and future reach of cosmic microwave background spectral distortion bounds. 
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I. INTRODUCTION 

Big bang nucleosynthesis (BBN) has been used for 
decades as a very powerful tool to constrain exotic parti¬ 
cle physics (for reviews, see for instance Refs. [H[3].) In 
particular, metastable particles populating the plasma in 
the early Universe could induce a nonthermal BBN phase 
via their decay products. Both hadronic and electromag¬ 
netic cascades typically contribute to these phenomena, 
although the former are much more model dependent. 
On the contrary, electromagnetic cascades are known 
to lead to a quasiuniversal q-spectrum, only dependent 
on the overall energy injected and the injection time: 
a monotonically decreasing, broken power-law (see e.g. 
Chapter VIII in Ref. [4] for a basic derivation). Recently, 
we pointed out that in a particular regime the commonly 
used universality of the spectrum achieved by photons as 
a result of electromagnetic cascades is violated pQ. This 
corresponds to the situation when the photons injected 
at energy E 1 are not sufficiently energetic to induce pairs 
onto the background photons at temperature T, and can 
be translated in the condition E 1 < 10 Ty ~ e y MeV (we use 
natural units with c = ks = 1). For T of order O(keV) 
down to (9(eV) characteristic of the period between the 
end of BBN to the formation of the cosmic microwave 
background (CMB), these energies are typically higher 
than the photodisintegration thresholds of light nuclei, 
denoted by E t h. As a result, the injection of photons of 
which the energy falls in the couple of decades from a few 
MeV to a few hundreds MeV may have an impact differ¬ 
ent than the one estimated with the universal spectra, 
given by 
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K 0 (f^) for E, < e x , 

Ko^y 2 for ex<E 7 <e c , (1) 

0 for E > e c . 


In the above expression, Ko = Eoe~£*[2 + ln(e c /ex)] -1 
is a normalization constant that enforces the condition 
that the total energy is equal to the injected electromag¬ 
netic energy, E 0 ; the characteristic energy e c = ra 2 /e™ ax 
denotes the effective threshold for pair production, e™ ax 
being the highest energy of the photon background onto 
which pairs can be effectively created; and ex e c /3 
is the maximum energy of up-scattered inverse Compton 
(IC) photons (See Refs. [5-7: for Monte Carlo studies 
leading to further justification of these parameters.) 

In our previous publication [T] , we illustrated how this 
may reopen the possibility of purely electromagnetic new 
physics solution to the so-called lithium problem, but we 
anticipated that other domains may be affected. Here we 
outline the impact on the constraints in the abundance 
vs lifetime plane for unstable early Universe relics, decay¬ 
ing electromagnetically, and derived from the deuterium, 
4 He and 3 He measurements. Our main conclusion is that 
the bounds are non universal and that they may be sig¬ 
nificantly more stringent than commonly thought. In the 
following, we will compare the constraints obtained from 
different elements in the hypothesis of the universal spec¬ 
trum with the actual constraints obtained for monochro¬ 
matic photon injections at different energies, below the 
pair production threshold e c . This parametrization is 
used solely for the sake of clarity; the differences would 
persist for any spectrum (either primary photons or sec¬ 
ondary due e.g. to upscattering of background photons 
via the IC by energetic e ± ) injected below the critical 
energy. 

This article is structured as follows. In Sec. [TT| we 
describe the features of the electromagnetic (e.m.) cas¬ 
cades and the breakdown of the universal nonthermal 
spectrum, as well as our method to solve the relevant 
Boltzmann equations. In Sec. |III[ we describe the non¬ 
thermal nucleosynthesis formalism and the observational 
constraints being used in the following. In Sec. [IV] we 
review the constraints coming from the CMB, notably 
from its spectral distortions, to which we will compare 
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the BBN ones. Our results are reported in Sec. [V] Fi¬ 
nally, Sec. [VT| contains a discussion with our conclusions. 


II. E.M. CASCADES AND BREAKDOWN OF 
UNIVERSAL NONTHERMAL SPECTRUM 

In general, in order to compute the nonthermal photon 
spectrum which can photodisintegrate nuclei, one has to 
follow the coupled equations of both photon and electron- 
positron populations. For the problem at hand, however, 
where we limit ourselves to inject photons incapable of 
pair production, it is a good first approximation to ig¬ 
nore the nonthermal electrons; while the injected pho¬ 
tons will in general Compton scatter and produce them, 
a further process, typically IC onto the photon back¬ 
ground, is needed to channel back part of their energy 
in the photon channel. The energy of these photons is 
significantly lower than the injected photon one: when¬ 
ever they are reinjected below nuclear photodissociation 
thresholds they are actually lost for nonthermal nucle¬ 
osynthesis; otherwise they would contribute to strength¬ 
ening the bounds, although only by a few percent, for the 
cases discussed below. For simplicity, let us also start by 
assuming that all photon interactions are destructive ; i.e. 
photons are not downscattered to a lower energy. Within 
this approximation, the Boltzmann equation describing 
the nonthermal photon distribution function / 7 reads 

EEA = -Y^,T{t))r{E^T{t))+S{E^i), ( 2 ) 

where S(E 1 , t) is the source injection term, T 7 is the total 
interaction rate, and we neglected the Hubble expansion 
rate since interaction rates are much faster and rapidly 
drive / 7 to a quasistatic equilibrium, = 0 . Thus, 

we simply have 
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where the term S for an exponentially decaying species 
with lifetime rx and density nx{t ), of which the total 
e.m. energy injected per particle is F’o, can be written as 
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and 0(100) MeV, and at temperatures of order few keV 
or lower, hence the thermal broadening is negligible, and 
a Dirac delta spectrum as the one above is appropriate. 

We calculate T 7 by summing the rates of processes that 
degrade the injection spectrum, namely: 

• Scattering off thermal background photons, jth- 
7 + 1th —^7 + 7- 

This has been studied in Ref. [ 8 ], and the scatter¬ 
ing rate of a 7 -ray with energy E^ over a blackbody 
with temperature T is given by 

3 


r 77 = —0.1513<a 4 ra e ( 



• Bethe-Heitler pair creation : 7 + N —>■ e ± + N. 
The cross section for this process is given by [9 
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• Compton scattering over a thermal electron : 7 + 
4 ^ 7 + e± - For the temperature of interest of 
O(keV), one can consider electrons to be at rest. 
In this case, we have [ 6 ] 
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where x = an( 4 n e is the number density of 
background electrons and positrons. 


In reality, not all scattered photons will be “lost”: even 
ignoring the energy transferred to e~ and e + , Compton 
scattering and 77 scattering still leave lower-energy pho¬ 
tons in the final state. This effect can be accounted for 
by replacing the rhs of Eq. ([2| by the following term: 

rOO 

S(E 7 , t) S(E 7 , t) + / dxK 1 (E 1 ,x,t)f 1 (x,t). ( 8 ) 

Je 7 


The additional term of which the kernel is K accounts 
for scattered photons and is obtained by summing the dif¬ 
ferential rates for the 77 scattering off background pho¬ 
tons and the Compton scattering over thermal electrons, 
respectively given by [ 8 ] 


with z(t) being the redshift at time t and the energy 
parameter (x (conventionally used in the literature) is 
simply defined in terms of the initial comoving density of 
the X particle n x and the actual one of the CMB, n 7 , 
via n x = n^(x/E q. We shall use as a reference spectrum 
the one for a two body decay X —)► 7 U leading to a 
single monochromatic line of energy Eq, corresponding 
to p 1 (E 1 ) = S(Ery — E 0 ). If the unspecified particle U 
is (quasi)massless, like a neutrino, one has Eq = mx/2, 
where mx is the mass of the decaying particle. Note that 
here, we will be interested in masses mx between a few 
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FIG. 1: Spectrum computed in this article (solid black line) 
compared with approximated one used in the literature (long- 
dashed red line), for the case Eq — 70MeV at T — 100eV. 
The short-dashed, blue lines show the contribution from the 
photon population as computed in our iterative treatment, 
with the number of iterations increasing from 1 to 7 from 
bottom to top. The dot-dashed curve is the estimated con¬ 
tribution to the photon spectrum from the nonthermal popu¬ 
lation of electrons excited by the energy loss channels of the 
photons. 


The integral in Eq. Q now depends on / 7 . We numer¬ 
ically solve this Boltzmann equation using an iterative 
method: we start from the Dirac distribution and the 
algebraic solution of Eqs. © and plug in the result 
thus obtained in Eq. © to estimate the new “effective” 
source term, and proceed. Note that the zeroth-order 
solution of Eqs. © and © is exact at the end point 
E 1 = Eq, with further iterations essentially improving 
the description at lower and lower energies. We stop 
iterating when the resulting improvement on the con¬ 
straints is smaller than 3%. Figure [l] shows the resulting 
spectrum proportional to / 7 according to the prefactor 
of Eq. Q] for an injected monochromatic photon of 70 
MeV at the temperature T = 100 eV in the commonly 
used universal spectrum approximation (long-dashed red 
line) and for the actual solution of the Boltzmann equa¬ 
tion, as a function of the iteration (short-dashed blue 
lines). For this example, one can estimate e c ~ 100MeV 
and ex — 30 MeV. Two features are clearly visible: i) 
the universal spectrum grossly fails for E 1 > ex, as ex¬ 
pected, since it imposes an artificial suppression; ii) the 
exact solution is significantly harder at intermediate ener¬ 
gies, but attains the same slope as the universal spectrum 
at low energies. However, the low-energy normalization 
is altered, since the universal spectrum is unphysical in 
pushing too many photons to low energies (below nuclear 
thresholds). 

Although neglected so far, an analogous treatment can 
be applied to the nonthermal electron distribution f e : 


the source term will be given by the Compton scatter¬ 
ing and Bethe-Heitler processes of nonthermal photons, 
and the “loss term” into the photon channel essentially 
by IC scattering. The latter will in turn correspond to 
a new source term in Eq. (J8|, the impact on the pho¬ 
ton spectrum is reported in Figure [l] with a dot-dashed 
green line. It is clear that, unless the injected photon en¬ 
ergies are too high, this only brings a modest correction 
to the low-energy tail of the spectrum, with the expected 
improvement in the constraints being even less promi¬ 
nent. The iterative solution technique adopted above 
would still perform correctly, although a detailed evalua¬ 
tion would render the calculation unnecessarily lengthy, 
and will not be pursued further here. We checked in 
fact that, for the cases discussed in the following, four 
iterations for the photon spectrum are enough to obtain 
bounds accurate at the 10% level (and often better) and 
always on the conservative side. 

Since the critical energy for pair production is a dy¬ 
namical quantity that increases at later times due to the 
cooling of the Universe, it may happen that the primary 
photons energy Eq is above threshold for pair production 
at early times and below it at late times (we do take into 
account that the decay is not instantaneous). In gen¬ 
eral, at each time we will compare Eq with e c and use 
the universal spectrum when Eq > e c or the monochro¬ 
matic spectrum with the complete expression for <S when 
Eq < e c . This gives always a qualitatively correct solu¬ 
tion, albeit it is somewhat approximate when Eq ~ e c . 
Since this is realized only in a very narrow interval of 
time, however, the final results are also quantitatively 
robust, barring artificial ”fine-tuned” results in a specific 
region of the parameter space. 


III. NONTHERMAL NUCLEOSYNTHESIS 
A. Review of the formalism 


At temperatures of few keV or lower, the standard 
BBN is over, and the additional nucleosynthesis can be 
simply dealt with as a postprocessing of the abundances 
computed in the standard scenario, for which we use the 
input values from Parthenope m , with the updated value 
of coming from Ref. m ■ 

As long as the amount of injected energy is small com¬ 
pared to the density of background CMB photons, one 
can neglect its impact on the expansion history. Thus, 
the nonthermal nucleosynthesis due to electromagnetic 
cascades can be described by a system of coupled differ¬ 
ential equations of the type 


dY A 

dt 
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where: Y A = TiA/^b is the ratio of the number density 
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of the nucleus A to the total baryon number density n 5 
(this factors out the trivial evolution due to the expan¬ 
sion of the Universe); and ct 7+ t^a is the photodissoci¬ 
ation cross section of the nucleus T into the nucleus A, 
i.e. the production channel for A; a 7+A ^p is the anal¬ 
ogous destruction channel. Both cross sections are ac¬ 
tually vanishing below the corresponding thresholds. In 
general one also needs to follow secondary reactions of 
the nuclear byproducts of the photodissociation, which 
can spallate on or fuse with background thermalized tar¬ 
get nuclei, but none of that is relevant for the problem at 
hand. According to Ref. m, the only signification sec¬ 
ondary production is that of 6 Li. Despite extensive work 
in the past, the current observational status of 6 Li as a 
reliable nuclide for cosmological constraints is doubtful, 
given than most claimed detections have not been ro¬ 
bustly confirmed, and a handful of cases are insufficient 
to start talking of a “cosmological” detection, see m • 
We shall thus conservatively ignore this nuclide and the 
secondary reactions in the following. 

With standard manipulations, namely by transforming 
Eq. (11) into redshift space, defining H(z) = H®( 1 + z) 2 
as appropriate for a Universe dominated by radiation 
(with H® = Hq Hq and being the present Hub¬ 

ble expansion rate and fractional radiation energy den¬ 
sity, respectively), one arrives at 

dY A = -1 

dz H®(z + l ) 3 

x ^ ^ Yj 1 I dE 7 f 7 (E 7 , >-a(-^7) 

- r J' dO 

poo 

- YaV / dE 1 f 1 (E 1 ,z)(J 1+A ^p(E 1 ) (12) 

P Jo 


which is solved numerically for the cases of interest. 


B. Light element abundances 

Among light elements, we can broadly speak of an 
agreement of standard BBN predictions with observa¬ 
tions for the case of 4 He, 3 He, and 2 H, while at face value 
the 7 Li yield is overpredicted by a factor ~ 3 with respect 
to observations. Since the interpretation of 7 Li observa¬ 
tions in terms of a primordial yield is still a subject of 
debate, see Refs. m m ns], one can consider two possi¬ 
bilities: either the observed values are not representative 
of the cosmological production mechanism, in which case 
it would be meaningless to derive constraints based on 
those observations, or alternatively, modifications to the 
standard BBN scenario, including electromagnetic cas¬ 
cades, could reconcile the envelope of 7 Li observed values 
with a primordial origin. In our previous paper [I], to 
which we address for further details, we discussed the 
latter possibility. In the following, we will adopt the for¬ 
mer, more conservative option, and hence we will not use 
7 Li for constraints on e.m. cascades. 


For the observationally imposed limits, we use the fol¬ 
lowing: for 4 He, which can only be destroyed by non- 
thermal BBN, we just impose the 2-cr lower limit on the 
mass fraction Y p > 0.2368 from Ref.[l4|. For 2 H we 
adopt the 2 -cr limit 2.56 x 10 -5 < 2 H/H < 3.48 x 10 -5 
from Ref. [15] : similar results would follow by adopt¬ 
ing the combination value compiled in Ref. [ 2 ], namely 
2.45 x 10 -5 < 2 H/H < 3.31 x 10 -5 , which is also closer 
to the results of Ref. PS; our interval also overlaps with 
the recent determination in Ref. El- For 3 He we im¬ 
pose no observational lower limit, but the 2 -cr upper limit 
from [E] 3 He/H <1.5x10 5 . It is worth noting that, 
had we used some alternative recent determinations such 
as m for 4 He or [20] for 2 H, some mild tension with the 
standard BBN predictions for the value of r] recently re¬ 
ported, e.g., by Planck would have appeared. These dis¬ 
crepancies are much smaller than the one affecting 7 Li, 
and could be easily accommodated with a more conser¬ 
vative error attribution: for 4 He this is the conclusion 
supported, e.g., in Ref. m or the recent Ref. m , essen¬ 
tially consistent with the value we quoted above; for 2 H 
it is also a possibility suggested by the slightly anomalous 
dispersion of several measurements around the mean (see 
e.g. the discussion in Ref. [ 2 ]). Alternative possibilities 
to reduce the tension with one or several of these deter¬ 
minations include a slightly different value of 77 between 
the BBN epoch and the CMB one, the addition of exotic 
phenomena such as cascades, and possibly others. In the 
following we shall adopt a similar attitude to the one 
adopted before for 7 Li and consider conservatively the 
more generous observational ranges reported above. This 
is also justified to ease the comparison with earlier anal¬ 
yses of cascade nucleosynthesis bounds, which used sim¬ 
ilar ranges. Our main emphasis here is in fact to gauge 
the impact of a more correct treatment of electromag¬ 
netic cascades, rather than deriving the most aggressive 
bounds achievable. Needless to say, should more precise 
observational values be confirmed in future studies, if in 
agreement with standard BBN expectations, it would be 
worth it to derive updated stringent bounds; if not, it 
would be interesting to rediscuss possible explanations 
in the context for instance of cascade nucleosynthesis, as 
we did for 7 Li in Ref. pQ. 


For the current application, the network of reactions 
used is reported in Fig. [2] and follows the parametriza- 
tion in the appendix of [ 12 ]. [Actually, the reaction 
4 He( 7 , 2 H) 2 H is significantly suppressed with respect to 
the others and thus not shown in the figure but is in¬ 
cluded in our numerical treatment.] Note that all cross 
section share the same qualitative features: they rise fast 
just above threshold, go through a peak (the so-called 
giant dipole resonance), eventually showing a decreasing 
tail at higher energies. We shall compare the bounds thus 
obtained with the ones coming from CMB spectral dis¬ 
tortions and entropy production, briefly recalled in the 
following section. 
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FIG. 2: Cross sections for the relevant photodisintegration 
processes. 


IV. CONSTRAINTS FROM THE CMB 


It is well known that a late injection of photons in the 
thermal bath can lead to additional measurable cosmo¬ 
logical alterations. 

For instance, the injection of a significant amount of 
energy can lead to modification of the photon-baryon ra¬ 
tio 77 or equivalently to the increase of the comoving en¬ 
tropy. Since the inferred values of U5 at the BBN and 
CMB epoch are compatible, no major entropy release 
could have taken place between nucleosynthesis and de¬ 
coupling. It can be shown that, in a Universe dominated 
by radiation and by considering that the decays have hap¬ 
pened at t ^ r, we have for a small fractional change in 
entropy (see e.g. Ref. m) 


^ =2.14x 10“ 4 

D Di 


Cx 




10- 9 GeV / V 10 6 s 


1/2 


(13) 


with a slight abuse of notation since (x^ now has to 
be intended to include any electromagnetically interact¬ 
ing decay products, all of which contribute to modify 
the photon-baryon ratio. To derive a statistically sound 
constraint, one should combine BBN and CMB data, al¬ 
lowing for an entropy increase between the two epochs. 
Since, as we shall see, this constraint is typically much 
weaker than others, such an exercise would bring us far 
beyond the scope of this paper adding a lengthy and un¬ 
necessary complication. We shall thus limit ourselves to 
illustrate the constraint that would follow by allowing 
a maximal 2% increase in the entropy between the two 
periods. This is an educated guess of the order of the 
bounds that one can expect, roughly corresponding to 
the 2-cr error bars on U5 from Planck 2015 m- 

Furthermore, as reviewed in detail in Ref. [ 22 ] . the 
spectrum of the CMB itself can also be affected through 
two types of deformation: a modification of the chemical 


potential /i and a modification of the Compton-7/ param¬ 
eter, which is related to the energy gained by a photon 
after a Compton scattering. To first order, it is possi¬ 
ble to distinguish the era of /i distortion from the era of 
Compton-// distortion, because the rate of the processes 
which are responsible for one type of distortions domi¬ 
nates at very different time. Basically, the [i distortion 
arises from rare process, implying a change in the num¬ 
ber of photons such asy + e—xy + y + e, whereas the 
Compton-7/ distortions are due to the end of the equilib¬ 
rium of Compton reactions, which happens much later 
in the history of the Universe, with a schematic /i — 7/ 
transition happening at for z ~ 4 x 10 5 , i.e. r ^5x 10 10 
s. For the relatively early time we focus on, the con¬ 
straints come essentially from /i -type distortions. We 
follow here the results of Ref. [ 22 ] . which contains im¬ 
provements with respect to the ones given in Ref. [24] , 
notably for 2 < 2 x 10 6 , while Ref. [24] is accurate enough 
at late times (see Fig. 16 in Ref. [22]). Hence, we adopt 


1/2 

11 ~ 8.01 x 10 2 ( ^ ) xjx 


1 GeV/ 


(14) 


with 


J = 


exp 




for z < 2 x 10 6 


2.082 ( ^ 


\ rx 


exp 


-1.988 


otherwise, 


(15) 


where r^c = 1-46 x 10 8 (To/2.7 K) 12 / 5 (Ot,/i 2 ) 4 / 5 ( 1 — 
V2) 4 / 5 is the “double Compton” interaction time in 
terms of the current CMB temperature To, with Y p ~ 
0.25 the primordial mass fraction of 4 He. We use the 
limit given by COBE on the chemical potential: |/x| < 
9 x 1CU 5 m, but we will also show the sensitivity that 
should characterize the future experiment PIXIE, of the 
order of fi > 5 x 10 — 8 , at 1-a [26] . 


V. RESULTS 


One of the most peculiar features of the spectral 
nonuniversality of photons injected below the pair pro¬ 
duction threshold is that the final outcome reflects the 
energy distribution of the injected photons with respect 
to the shape of the relevant photodisintegration cross sec¬ 
tions, shown in Fig. [ 2 ] This motivated us to choose in 
the following for each nuclide, the results for two repre¬ 
sentative examples of monochromatic injection: one close 
to the resonant peak and another one well after it. The 
markedly different outcomes obtained in the two cases 
should thus convincingly argue that constraints of actual 
models are going to be determined not only by the de¬ 
cay time and the overall energy injected but also by the 
energy range at which the bulk of the photons lies. 
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A. Constraints from 4 He 

The simplest situation is certainly the one concerning 
4 He: being the only abundant nucleus subject to photo¬ 
disintegration, its nonthermal e.m. production is irrele¬ 
vant, and one only has to care about its destruction; i.e. 
only the term proportional to Ya at the rhs of Eq. © 
is important. The results obtained by using a monochro¬ 
matic injection at 70 MeV (hatched/light shaded red), 
at 30 MeV (dark shaded red), and the universal spec¬ 
trum are shown in Fig. [3] The vertical lines indicate the 
time at which the threshold energy for pair production e c 
starts exceeding the corresponding injected energy. One 
might naively expect that this is the time at which the 
constraints obtained from the incorrect use of the uni¬ 
versal spectrum start to deviate from the actual ones. 
However, when taking into account the fact that the de¬ 
cay is not instantaneous, it turns out that constraints 
already start to deviate at ~ Tx / 5, and the closer to the 
post-threshold cross section resonance we inject energy, 
the earlier deviations appear. 



FIG. 3: Constraints from 4 He depletion in the standard 

case (black line) and for a non-universal spectrum with Eq — 
30 MeV (dark shaded red) and Fo=70 MeV (hatched/light 
shaded red). We also show the sensitivity to the entropy vari¬ 
ation constraint (green dashed line), current constraints from 
CMB spectral distortions (excluded above the short-dashed 
blue line), and the sensitivity reach of the future mission 
PIXIE 26] (above the red dot-dashed line). 

Although the BBN bounds coming from excessive de¬ 
pletion are typically (but not always!) weaker than CMB 
distortion bounds, also reported in the figure, note that in 
both cases bounds can differ from the ones derived with 
the universal spectrum by a large factor, up to an order 
of magnitude if the energy of the photons is around the 
peak of the photodissociation cross section. For higher- 
injected energies, they tend to become closer to the uni¬ 
versal spectrum constraints, as it should. In fact one 
can envisage fine-tuned situations in which they become 



FIG. 4: Constraints from deuterium depletion in the stan¬ 
dard case, with Eq < 20 MeV (black solid line) and for a 
nonuniversal spectrum with Eq — 4 MeV (dark shaded red) 
and Eq—20 MeV (hatched/light shaded red). All other con¬ 
straints/sensitivities shown as in Fig. [ 3 ] 

slightly weaker, albeit this conclusion does depend on the 
extrapolation of the photodisintegration cross sections, of 
which the reliability at high energy has never been quan¬ 
titatively assessed in the context of BBN applications. 

B. Constraints from 2 H 

For deuterium, the situation is more complicated be¬ 
cause several regimes are present. At low tx, e c is below 
4 He photodissociation threshold (and in some cases also 
below A = 3 nuclei photodissociation thresholds, which 
are, however, less relevant). Hence, only constraints from 
overdestructions are present. At high tx, however, what 
dominates is the overproduction from 4 He destruction. 

Figure [4] shows the illustrative case where production 
channels are turned off: this is exact for Eq < 8 MeV, 
but a good approximation till Eq < 20 MeV. Note the 
qualitative similarity to the 4 He case, apart for the mod¬ 
ifications due to the different features of the respective 
cross sections. 

Whenever production channels from 4 He are open, 
which requires Eq > 20 MeV, the constraints are signif¬ 
icantly stronger at large tx, as shown in Fig. [5] Once 
again, a violation of universality (and a sensitivity to 
the energy dependence of the cross section) is clearly 
manifest by the two cases shown, Eq = 30 MeV and 
Eq = 70 MeV. 

It is also worth noting that for deuterium the con¬ 
straints are significantly stronger than the ones coming 
from CMB spectral distortions. By improving the sen¬ 
sitivity to /i down to /i ^ 5 x 10 -8 , the sensitivity ex¬ 
pected by the future mission PIXIE [26] (shown by the 
red, dot-dashed curve) would greatly strengthen these 
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FIG. 5: Constraints from Deuterium depletion and produc¬ 
tion in the standard case, with Eq > 20 MeV (black solid line) 
and for a nonuniversal spectrum with Eq — 30 MeV (dark 
shaded red) and Eo—7t) MeV (hatched/light shaded red). All 
other constraints/sensitivities shown as in Fig. [ 3 ] 



FIG. 6: Constraints from the 3 He production in the standard 
case (black line) and for a nonuniversal spectrum with E 0 — 
30 MeV (dark shaded red) and F?o=70 MeV (hatched/light 
shaded red). All other constraints/sensitivities shown as in 
Fig.0 

constraints, with the exception of very small lifetimes 
where deuterium overdestruction would still provide the 
dominant bounds. 


C. Constraints from 3 He 

First of all, a premise is necessary: there are in fact 
two nuclei with A = 3, 3 He, and 3 H, the latter being 
unstable to beta decay into 3 He with a half-time of over 



FIG. 7: Global BBN best constraints in the standard case 
(black line) and for a nonuniversal spectrum with Eo — 
30 MeV (dark shaded red) and F?o=70 MeV (hatched/light 
shaded red). All other constraints/sensitivities shown as in 

Fig. [| 

12 years, or about 4 x 10 8 s. Practically, however, for the 
purposes of the constraints discussed here, one can sum 
the equations for 3 He and 3 H and treat them as a sin¬ 
gle effective nucleus with A = 3. The reason is twofold: 
first, we only require 3 He not to be overproduced with 
respect to the observational upper limit. Hence, the key 
reactions are the production channels by single nucleon 
photodisintegration from 4 He, which are only open above 
20 MeV, rather than the destruction ones. Second, 3 He 
and 3 H are “mirror nuclei” under the isospin symmetry 
nop, and their nuclear properties are in fact very sim¬ 
ilar: the corresponding thresholds in nuclear cross sec¬ 
tions, for instance, only differ by some 0.8 MeV (compare 
the two curves in Fig. [2j) From Fig. [6j where we report 
our results, it is clear that the photodisintegration cross 
section for single nucleon emission from 4 He, when open, 
is so important that very stringent nucleosynthesis con¬ 
straints follow. In fact, they are much stronger than the 
current ones coming from CMB spectral distortions, al¬ 
though future PIXIE sensitivity might improve over them 
over most of the parameter space. 

Notice the importance of the nonuniversality: the two 
cases with 30 or 70 MeV monochromatic injections lead 
to significantly different constraints, and in both cases 
depart from the “universal spectrum” ones. 


VI. CONCLUSIONS 

We have argued that the universality of the photon 
spectral shape in electromagnetic cascades has often been 
used in cosmology even beyond its regime of applicabil¬ 
ity. When the energy of the injected photons falls be¬ 
low the pair production threshold, i.e. approximately 








































when E 1 < ra 1 2 /(22T) ~ lOT^yMeV, the universal 
form breaks down. In Ref. pQ, we showed how this could 
potentially open the possibility of a purely electromag¬ 
netic decay solution to the so-called lithium problem. In 
this article, we showed how important the modifications 
to the photon spectrum in this regime are for the con¬ 
straints from nonthermal BBN. This required the numer¬ 
ical solution of the relevant Boltzmann equations, which 
we attacked by an iterative scheme. 

The constraints we obtained, for illustrative cases of 
monochromatic energy injection at different epochs, are 
often much stronger than the ones presented in the lit¬ 
erature (up to an order of magnitude), notably when 
the injected photon energy falls close to the peak of the 
photodisintegration cross section of the relevant nucleus. 
In fact, the breaking of the nonuniversality is nontriv¬ 
ial and is essentially controlled by the energy behavior 
of the cross sections: in the universal limit, most of the 
photons he at relatively low-energies, so that the cross- 
section behaviour at the resonance just above threshold 
is what matters the most. In the actual treatment, the 
photons may be also sensitive to the high-energy tail of 
the process. Future studies aiming at assessing the nu¬ 
clear physics uncertainties affecting these types of bounds 
would benefit from this insight. It cannot be excluded 
that in some cases constraints weaken a bit with respect 
to what is considered in the literature. 

We also compared BBN bounds with constraints com¬ 
ing from CMB spectral distortions. A summary plot 
of the “best constraints” is reported in Fig. [7| for two 
choices of the monochromatic photon energy. We con¬ 
cluded that BBN limits are improving over current con¬ 
straints from COBE via the requirement not to under¬ 
produce 2 H (at low injection lifetime tx), or not to over¬ 
produce 3 He (at high tx), while 4 He is never competitive. 
The bounds from a future CMB spectral probe such as 
PIXIE would not only greatly improve current CMB con¬ 


straints but would also reach the level of current con¬ 
straints from 3 He (often improving over them) allowing 
for an independent consistency check. This is reassur¬ 
ing, since the cosmological reliability of 3 He constraints 
does stand on some astrophysical assumptions. Below 
tx ~ 5 x 10 5 6 s, however, 2 H constraints would probably 
remain the most stringent ones for a long time to come. 
Fortunately they are i) quite robust, relying on the single, 
well-known cross section 2 H(y,n)p, and ii) easy to com¬ 
pute, since no coupled network of equations needs to be 
solved, the problem reducing to the numerical evaluation 
of a single integral (the same situation leading to Eq. (6) 
in Ref. pQ.) This is also the region where constraints 
coming from hadronic decay modes (not revisited here) 
are quite stringent. A synergy between BBN and CMB 
is thus going to be necessary for this kind of physics even 
in the decades to come. 

In conclusion, our work suggests that models in the 
literature that fulfilled the BBN constraints with less 
than an order of magnitude margin should perhaps be 
reconsidered. In particular, those characterized by soft 
gamma-ray emissions and/or at relatively late times 
should have been more prone to incorrect conclusions 
about their viability. Our study also suggests that ac¬ 
tual bounds should be derived via a case-by-case analy¬ 
sis. Finally, we provided further arguments supporting 
the usefulness of an improved constraint from CMB spec¬ 
tral distortion of the /a type, since it would not manifest 
the unexpected sensitivity to the shape of the energy in¬ 
jection that we have uncovered. 
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